An extended massless phase and the Haldane phase in a spin-1 
isotropic antiferromagnetic chain 



Chigak Itoi and Masa-Hide Kato 
Department of Physics, College of Science and Technology, Nihon University. 
Kanda Surugadai, Chiyoda-ku, Tokyo 101, Japan. 
(February 6, 2008) 



Abstract 

We study the phase diagram of isotropic spin-1 models in the vicinity of the 
Uimin-Lai-Sutherland (ULS) model. This is done with the help of a level-one 
SC/(3) Wess-Zumino-Witten model with certain marginal perturbations. We 
find that the renormalization group flow has infrared stable and unstable 

trajectories divided by a critical line on which the ULS model is located. 
The infrared unstable trajectory produced by a marginally relevant pertur- 
bation generates an exponential mass gap for the Haldane phase, and thus 
the universality class of the transition from the massless phase to the Haldane 
phase at ULS point is identified with the Berezinskii-Kosterlitz-Thouless type. 
Our results support recent numerical studies by Fath and Solyom. In the 
massless phase, we calculate logarithmic finite-size corrections of the energy 
for the St7(^)-symmetric and asymmetric models. 

PACS numbers: 75.10.Jm, 75.40. Cx, 05.70.Jk, ll.10.Hi, 11.25.Hf 
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I. INTRODUCTION 



The phase diagram of isotropic spin-1 chains has not yet been understood sufficiently. 
The characteristics of ground states can change drastically depending on a coupling constant 
of the model jl] . Even though there are many rigorous |2|j3[] and exact 0-0] works at several 
isolated points, one encounters unconformable issues in a certain region, especially in the 
non-integrable region around an integrable point of the Uimin-Lai-Sutherland (ULS) model. 

The general form of the spin-1 Hamiltonian which consists of nearest neighbor interac- 
tions with rotational symmetry is 



H{d) = £ \cos9 (5 j ■ S j+1 ) + sm6 (Sj ■ S j+1 ) 



(1-1) 



where the coupling constant is controlled by one parameter 9 G [0,27r). It is our main 
concern to understand the macroscopic behavior in the vicinity of the ULS point 9 = tt/A 
||||. It is known that the ULS model has massless excitations described by the Wess- 
Zumino-Witten (WZW) model. The region \9\ < 7r/4, which contains the standard Heiseberg 
antiferrromagnet (9 = 0), is believed to be in the Haldane phase which has only massive 
excitations, as suggested by some numerical works |8||| and rigorous studies at 

9 = tam 1 (1/3) [0,0. On the other hand, the natures of the model in the region tt/2 > 
9 > it / 4 are theoretically less understood. 

In this paper employing a renormalization group (RG) method in a continuum field 
theory, we show that the region 9 > 7r/4 near 9 = tt/A is a massless phase, and that the 
phase transition from this massless phase to the Haldane phase at ULS point belongs to 
the Berezinskh-Kosterlitz-Thouless (BKT) type universality class. This result is consistent 
with a numerical study obtained by Fath and Solyom jlgR . For this purpose, we map the 



ULS model to the SZ7(3)i WZW model, which reproduces some exact results obtained from 
the Bethe ansatz 110,0-0. In the non-integrable region around the ULS point, we show 
that the SZ7(3)i WZW model is perturbed by adding a S£/(3)-breaking marginal operator 
which causes the BKT transition. We observe several nontrivial behaviors as in some other 
CFT deformed by marginal operators [TP|--f2!2j| . Despite a number of studies on the BKT 



transition and the logarithmic corrections in SZ7(2) systems, those concerned with SU (u) 
symmetry for v > 2 has been seldom discussed. Here, we study the BKT transition and the 
logarithmic correction in the v > 2 case and we find its different nature from u = 2 case. 
The obtained continuum theory enables us to calculate the logarithmic finite-size correction 
in the energy of the ground state and the first excited states in the region of the massless 
phase. Following Ludwig and Cardy [p4] , p5|1 , the finite size correction to the ground state 



energy of the model in a strip space with the width L is 

1 \ (1-2) 



TTV 

^G.S = £ oo L ~ ^ C ( L )' 



m d G.s Ln fin (In L) 

c{L) = c vir + - + O 



(InL) 3 V (InL) 4 ' (lnL) 4 / ' 

where is the non-universal bulk contribution to the ground state energy depending on 
cut-off scale. The minimal energy of an excited state related to a certain primary field with 
conformal weight x n /2 is given by 
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£ n = £ G .S + -j^lni L ), 

. _ . oL ^ An (In L) 1 

7.(i) = *.+ taT + ^TCip^ 

where d n is a coefficient of a certain three point function. We calculate these universal 
coefficients of the logarithmic corrections by the obtained continuum field theory. 

The outline of this paper is as follows. In Sec. |TJ, a strong coupling abelian gauge theory 
is introduced as a critical field theory of the ULS model, which allows us to evaluate exact 
values of universal quantities. We show the equivalence of this critical theory to the level-one 
SU (3) WZW model. This argument can be generalized to a certain SU(v) symmetric spin 
model which includes the ULS model in the v = 3 case. In Sec. [TT1], we discuss an extended 
non-integrable spin model with S[/(^)-asymmetric interaction on the basis of the level-one 
SU (u) WZW model with an asymmetric perturbation. We pin down the marginal operator 
J2a=i ^LafiK 2 ) 3ro.b (X) m the SU{v)\ WZW model as the S[/(z/)-asymmetric interaction in 
the original spin model. 

The logarithmic corrections of its energy in the massless phase are evaluated and the 
difference between S?7(3)-symmetric and asymmetric model is indicated. Finaly, we discuss 
the universality class of the transition from the massless phase to the Haldane phase which 
belongs to the BKT type. 




II. CRITICAL THEORY OF SU(i/) SPIN CHAIN 



To begin with, we extend a QED 2 description for S£/(2) spin model [^J to the SU{y) 
one. The spin chain is mapped to the WZW model with some perturbations by this method. 
We redefine the Hamiltonian (|1 . 1|) near 9 = 7r/4 as 



cos 9 



H{o) = T,\{s r s, J+1 ) +1 (s r s 3+l 



3=1 



(2.1) 



with 7 = tan#. We use fermion operators cj a , cj a for the spin variables 



(2.2) 



a,/3=l 



where L x , L y and L z are spin-1 matrices. In this case, eq (|2.1|) can be expressed in the 
fermions [15] 



3=1 



(2.3) 



in which a trivial constant is neglected. Here the local constraint, X^ =1 Cj a c j,a = 1, is im- 
posed in order to restrict the dimension of physical space to three at each lattice site. Due 
to this constraint, empty, double and triple occupancy states of the fermions are forbidden 
at each lattice site. The first term in eq( |2.3|) is a exchange operator between nearest neibour 



3 



sites and the second one is a projector onto a singlet bond. Eq(|2~3|) has the local U(l) sym- 
metry, a translational symmetry by one lattice site for all values of 7. In addition, a global 
SZ7 (3) symmetry appears at the ULS point (7 = 1) where the Hamiltonian (j2.3|) consists 
of only exchange operators, and the model becomes Bethe ansatz solvable. This fermion 
expression can be extended to Bethe ansatz solvable model with higher spin. When spin-S", 
2S + 1 kinds of fermion on each lattice site are introduced, and the constraint on each site 
is given by J^^Ji 1 c },a c j,a = 1- In general, bond interactions of an isotropic spin-5 1 chain are 
represented 

by a polynomial of X = Si-Sj. These integral families for higher spin chains are classified 
by Batchelor, Yung and Kennedy |H5 |. One of those is the ULS model with an arbitrary 



spin. We give the expressions of the exchange operator in terms of spin matrix for S < 2 in 
Table |. Hereafter we discuss the fermionized Hamiltonian ( |2.3| ) as v species [v = 2S + 1). 

The euclidean action is written by introducing Lagrangian multiplier x f° r the local 
constraint and the Hubbard-Stratonovich transformation 

r/3 L 

A i = / dr E { c L d rCj, Q + iXj(c},« c i.a - !) + ?W(7)) , (2-4) 

3=1 

where (3 is an inverse temperature. The Hamiltonian for one bond interaction is 

^■3,3+l(l) = ~ C },aQj,j+l C j+l,a ~ C ]+l,aQj,j+l C j,a + (7 — ^) C ],a C j,/3 C ]+l,a C j+l,/3- 

The complex auxiliary fields {Qjj,}, {Qjj,} are introduced to decompose the two-body 
fermion interaction into single body. A local U(l) gauge transformation 

c jt a -> ' '"' ' ( 'j.r, ■ Xj -> Xj ~ d T (pj, Q,,< -> Q.j.j'> ! -' ; , (2.5) 



preserves eq([ 

First, we study SU(u) -symmetric point 7 = 1 by the mean field theory without taking 
into account the local constraint YZt=i c j,a c j,a = 1- We shall treat the local constraint later. 
In the mean field theory, the auxiliary field Qjj> is a constant Ro, and then dispersion 
relation becomes e(k) = — R Q coska, where a is a lattice spacing. The ground state is given 
by the fermi sea filled up to fermi level ±kp with hp = ir/ua. The low energy physics can 
be described in terms of ip L and tp R which is the lattice fermion operator Cj, a only around 
the fermi surface with a certain low energy cutoff A(<< hp) as 

—r c i^ — Aa( x ) ex P (~ik F x) + ip Ra (x) exp (ik F x), x = j a. (2.6) 



In this representation, the local gauge transformation eq( |2.5|) corresponds to the U(l) vector 
transformation 

A translation by one-site on the original lattice space, 

Cj, a -> c j+ltCe exp(ik F a), 
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corresponds to a chiral Z v transformation: 

(Aa, Vw) ex P (2iA*a)) . 

As far as the translational symmetry is not broken, the effective field theory becomes chiral 
.Z^-invariant. 

Now, we take into account the deviation from the mean field approximation. In the 
following parametrization of the auxiliary field 

Q jtjl = R(iUL)e Xp {i\j-j'\A 1 {^)}, 
the deviation of Qjji becomes 

Qj,j' ^R + SR(^) + ili - f\M ljj f)- (2-7) 
The local constraint is expressed as 

4aAM = Q, 4aAa(x)=0, (2.8a) 

C^W + C^W=0- (2.8b) 
We obtain a chiral Z'^-invariant effective Lagrangian in terms of low energy variables: 

C = C + C ini , 

where 

C = 2$ id + L4) i, La + 2^ Ra (d + L4) *P Ra , ^ 

C mt = constant x V^A/A^W 

Here, the gauge field A is a low energy variable corresponding to the Lagrangian multiplier 
X = aA and A = A + L4i,A = A — L4i. This effective theory is a perturbed abelian 
gauge field theory with a sound velocity v = R^a sin (kpa), here t> is set to unity. In deriving 
the Lagrangian L we have picked up the terms to 0(a 2 ) and neglected the highly oscillating 
terms and higher derivative terms. The four-fermi interaction C int is induced by performing 
the gaussian integration over the 5R - field and also by the second constraint ( |2.8b| ). This 
interaction can be expressed in the form 

-^)m+2 u j:j L A (z)j R A (z), 

V A=l 

where j L(R) = i l l(R )a i ! L{R)a and J^ R) = tpL(K)aT A ^ L{a) p- 

In Appendix A, the SU(u) basis is summarized. We should define a regularization for 
the U(l) current which preserves the local Uy(l) gauge symmetry ipL — > exp (ia)ip L , ipR — > 
exp — (ia)^. An arbitrary local composite operator should be defined in the gauge invariant 



point spritting regularization [27]. The current operators are defined in this way as well. 
We shall discuss the importance of the marginal perturbation £ int on the basis of a RG 
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calculation later. To solve this system £ ? we take a gauge fixing condition A = and a 
parametrization of A with a scalar field <j>(z, z) 

A = d(f)(z,z), ^(z)=^(*)exp(-i0). (2.10) 

The gauge invariant regularization defines the unique Jacobian for the chiral transformation 
ip L — > ip L which induces the kinetic term of the scalar field . Then, the two-dimensional 
abelian gauge theory with global symmetry SU(u) is expressed as a decoupled free field 
Lagrangians: 

f d^z 

/ T> (^matter •^-'gaugo ^--ghost) (2'H) 
J Z7T 

with 

^matter = ^IcfiAa + 2 V4^VW> 

^gaugo = -vdcfdcj), (2.12) 
-Cghost = 2r/<9e + 2r]de, 

where the last term is the Fadeev-Popov ghost one originating from the measure of the 
gauge field parametrizedin terms of </>. Hereafter, the tilde of the left moving fermions will 
be omitted for simplicity. The operator product expansions (OPE) between free fields are 

Z — Uu 
Z — uu 

e i<t>(z,z) e -i<t>(u,u>) „ | z _ + . . . ^ 

which allow us to calculate OPE for energy-momentum tensors and read off the central 
charges 

Cmatter Cgauge 1, Cgh os t 2. 

As expected, the total central charge c to tai equals to l v — V which agree with the Bethe 
ansatz's result fl7fl . The negative sign in the Lagrangian £ ga uge, eq (|2.12|) , suggests that the 



U(l) degrees of freedom freeze in the asymptotic behaviors of the spin system. Since the 
conformal weight of the vertex operator becomes negative, eq( |2.13|) shows unphysical infrared 
behavior, and thus it should not appear by itself. Actually, the U(l) current regularized in 
the gauge invariant way 

(2-14) 

( |2.14| ) has no Goto-Imamura-Schwinger term in the U(l) Kac-Moody algebra. Therefore, 
this U(l) Kac-Moody algebra has only a trivial representation j L = 0, j R = 0. According to 
the bosonization formula in the U(z/) = U(l) x S?7(^)-invariant free Dirac theory, degrees 
of freedom of U(l) (charge) and SU(u) (spin) are separated into a scalar boson and a SU^u) 
WZW theory. In our case, the U(l) degree of freedom is killed by the gauge field which is 
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represented by the scalar boson with negative norm. Eq (|2.11|) can be identified with the 
level-one SU(u) WZW model [28]. The Wess-Zumino primary field G(z,z) is given by 

G ap (z,z) oc 4W^%(2), G(z,z) G SU(u) 

and the conformal weight is {y — l)/2u. To compute the asymptotic behavior of the spin 
correlation function for the bulk, it is enough to replace the spin-one operators ( |2.2p by the 
continuum fields in terms of eq( |2.6| ) and eq (|2.10| ) 



a- 1 S i ~J L (r) + J R (r) + 



where J L(R) {r) = ^ L(R)a (r)L a ^ L(R)l3 
of Tomonaga-Luttinger liquids |29[ 



^ (r) iWi^(r) exp (2ik F r) + h.c , 
r). Using this, we obtain the typical correlation function 



. _ _ . 1 cos (2k F r) 

{S r ■ b ) oc — - + constant x - 



.2.1- 



(2.15) 



with scaling dimension x = 1 — 1/u, in which the second term is dominant as r — > oo 
and the momentum distribution shows a power law singularity near the fermi momentum 
hp. The appearance of the oscillating factor is reflection of the chiral Z v symmetry in the 
ant if er romagnet . 



III. THE ROLE OF MARGINAL OPERATORS 

We have neglected the marginal operators so far. One of them is the SU(u) current 
interaction in C mt which gives logarithmic finite-size corrections. Besides, there is another 
operator J2j=i c ],a c j,/3 c ]+i,a c j+i,/3( x ) which breaks the global SU(u) symmetry except at the 
ULS point. The continuum form of the S[/(^)-asymmetric interaction is given in terms of 
eq(p|) by 

C j,a C j,f3 C j+l,a c j+l,f3 — ^LaVh^Ra^R^ + ' ' ' • 

which is also chiral Z^-invariant. The corresponding field theory is expressed by the WZW 
model with these marginal perturbations without global SU(v) symmetry. This SU(u)- 
breaking operator becomes marginally relevant for the coupling constant 7 < 1, and thus a 
dynamical mass generation is expected. 

We consider a perturbed CFT with the following action: 

A = A su{u)i +j:g t J^H^), (3-1) 

where 

yv 1 — 1 
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There are no other relevant or marginal operators with rotational and chiral Z% symmetry. 
The coupling constant g 2 is proportional to 7 — 1 with a positive coefficient in the case 
of v = 3 (S = 1). The unperturbed action A su , v \ is given by eq( |2.11| ) and the marginal 
operators obey the OPE algebra 

^ l \z,z)^(o,o) ~ -L _ _L$(i)( ,o) + ■ • • , 

\z\ \^\ 

^ (z, z)^ (0, 0) ~ -L + -^$( 2 ) (0, 0) + ■ ■ • , 

Z)*< 2 >(0, 0) ~ - -Aj ($ (1) (0, 0) - *< 2 >(0, 0)) + ■ • ■ , 

(z/ + 1) \z\* \ z \ 

where 

2z/ r 2 



This algebra gives the following one- loop (3 functions 

$2(91,92) = = ~92 2 -hi92 + (92 3 i9i 2 ), 



(3-3) 



where e' = a. These coupled differential equations can be solved in an integral form thanks 
to a conservation law. An arbitrary trajectory in the coupling constant space (gi, g 2 ) obeys 
the following equation: 

X 2 -Y 2 = C\Y\^ v - i)tv (3.4) 

with X = g 1 — g 2 and Y = —g 1 — g 2 , where C is an arbitrary real constant. The sign 
of the initial value of g x should be chosen to be negative in order to agree with the result 
of the Bethe ansatz. The running coupling constant g 1 is renormalized to be zero in the 
infrared limit. Therefore this model is in the second region g x < 0,g 2 > 0or the third one 
g x < 0,g 2 < in the coupling constant space. In the v — 2 case, this perturbed CFT 
describes the well-known spin-1/2 XXZ chain and eq( |3.4j ) shows a hyperbolic trajectory 
where the BKT transition occurs beyond the SC/(2) symmetric line X ± Y = 0. There is 
one parameter family of fixed points in c vir = 1 , that is a fixed line g\ + g 2 = [33|. Note 



that the topology of the flow diagram in the case of v 7^ 2 differs from that in v = 2. The 
only fixed point is = g 2 = except for v = 2. 

Here, we show main results of the RG flow which will be illustrated in the remaining 
part of this paper. The RG argument classifies the coupling constant space with g± <0 into 
the following three cases: 

• #2 = 0; SU(u) symmetric and asymptotically non-free, 

• g 2 > 0; SU(u) asymmetric and asymptotically non-free, 

• g2 < 0; SU(u) asymmetric and asymptotically free. 
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Since the interaction <3>( 2 ) is marginally relevant for g 2 < 0, and is marginally irrelevant for 
#2 > 0, the trajectory along g 2 = becomes the BKT transition line. In a finite system in the 
asymptotically non-free region g 2 > 0, thermodynamic quantities acquire some corrections 
due to the presence of marginally irrelevant operators, while in an infinite volume limit there 
is no influence from them. We indicate the difference of the finite-size corrections between 
SU(v)- symmetric and asymmetric models. In the third case of g 2 < 0, the marginally 
relevant interaction $^ 2 - ) can generate a mass gap which might be interpreted as the Haldane 
gap. 



First, following Ludwig and Cardy [24,25], we calculate finite-size corrections in the 



SU(v) symmetric model (g 2 = 0) with g\ < 0. The finite-size corrections to the ground state 
energy of the S?7(z/)-symmetric models in the eq( |1.2| ) are calculated as 

c vir = v-l, d GS = (3.5) 

where v = 2S + 1 for the spin-S*. 

The finite-size corrections to the low-lying excited energies are calculated from the most 
relevant primary field 

O a (z,z) = 4 a (z)T^(z)e^ +h.c, 

where T A 's for A = 1, ■ ■ ■ ,u 2 — 1 are the SU(v) basis, T° = I/^/2u, and they are also 
normalized as Tr[T A T B ] = S AB /2. The primary states, \O a ) = \im z ^ O A (z,z)\0) , become 
eigenstates of Virasoro's charge L (L ) with an eigenvalue x/2. Their OPE are given by 

O A (z,z)O B (0,0)~^ z ^ + ---, 



\z 



O A (z,z)$ (1) (0,0) ^0^(0,0) + 



\z 



with the OPE coefficients 

1 f u 2 - 1 for A = 0, 

= , X < 

v\/v 2 - 1 | -1 for A = l,...,z/ 2 - 1. 
We obtain the universal quantities in eq (|L3| ) 

26, f 1 - \lv 2 for A = 

, x = i-i/„, rf . = ^ = { _ iy imA = h ... y _ h (3-6) 

These v 2 states are classified by the total spin. As shown in Appendix A, the state A = 
describes the singlet excitation and other v 2 — 1 primary states are higher spin states with 
spin up to (v — l)/2. In the finite-size corrections up to the logarithmic size dependence, 
the singlet excitation is not favored compared to those with higher spin. 

The finite-size effect for the spin correlation function (|2.15| ) is obtained immediately from 
the information on the excited energy [pi|, P^ , P^ ] . 
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where 



(5 r • 5 > « 008(2^)^(^(0,^, G A (g x (r),r) = -S^, (3.7) 

7" A 



2 

(74 = -2d A = — 



except for A=0. Our results for v = 2 listed in Table |T| agrees with the Bethe ansatz's ones 
in ref [^,^]. The leading finite-size corrections c V i r and xa in the SU(^) symmetric model 



agree with Bethe ansatz [O,?], as well. 



Now we consider the second case g 2 > 0, where there is a marginally irrelevant 817(1/)- 
asymmetric interaction. The situation is crucial whether v = 2 or not. Even though the 
action describing the ultraviolet theory has no SU(u) symmetry due to the SZ7 (zz)-breaking 
interaction, the SU (v) -breaking interaction have no effect on the leading terms of the finite- 
size correction except in v = 2 case. The RG indicates that the SU(u) symmetry appears 
dynamically for macroscopic scale even though the g 2 -term in eq( |3.2| ) is switched into the 
fixed point action. The difference between the S£/(z/)-symmetric and asymmetric model 
appears in the logarithmic correction term. 

To calculate the logarithmic correction, we note that the RG flow eq( |3.3|) with an initial 
condition g 1 < and g 2 > is absorbed into the fixed point along the line 

9i = ~92- (3.8) 



The macroscopic property of the system is determined by the scale I ^> 1, and we can 
estimate the deviation from the line as \gi(l) +^(0 1 ~ 0(/ _ ~) with the help of the integral 
curve eq ( |3.4|) for arbitrary solution with an initial condition in the second region. Therefore 
we can calculate the logarithmic correction by assuming that the marginally irrelevant flow 
for g 2 > is described by the action 

d 2 z 



A = A su[u)i + gi J 2^*0,2) 



with \l/(z, z) = \J~^f (^\z, z) — & 2 \z, z)\ which is normalized by 



\zr \zr 



where the OPE coefficient B is 



B=1 fc±i(6-26). 



2v 

This assumption might hold, since the current of the RG would spend the fair time near the 
fixed point with dilatation. As in the discussions of the symmetric model, we can evaluate 
the coefficients of the finite-size energy correction from the one-loop renormalization which 
obeys dg x jdl = (B/2)g 1 2 . For the ground state energy, we obtain 

v(v-l) 

C vir — V ~ 1) "G.S — 77J _ 2^2' 
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in which the logarithmic coefficient is different from eq( |3.5| ). The three-point function in the 
expression of the excited energy is given by using the OPE 



O A (z,z)®(0,0) ~ -^-O A (0,0) + 

\zr 



Here the coefficient B A takes three different values according to the symmetric properties of 
the matrices {T A } under the matrix transposition. These are given by 

( v - 1 for A = 
B A = ^=^= x J 1 for 0) with \T A ) = — T A 
yj2v(v-l) [ _i for A(y^ 0) with (T ) = T A . 

As a result, we have the universal coefficients in the anomalous dimension (11.31) 



v - 1 for A = 

r,: i-l/i/, d A = ^^x<{ 1 for 0) with *(T A ) = - T A (3.9) 



i/-2 



-1 for 0) with t (T j4 ) = T A . 



The OPE coefncents B and _B A give the exponents <J A ^ characterizing logarithmic distance 
dependence in eq (|3.7| ) 

The primary states with A ^ in the symmetric model are degenerate even if we consider 
the logarithmic correction, those in the asymmetric one splits to two levels. As shown in 
Appendix A, the difference of the OPE coefficients because of the symmetric and antisym- 
metric properties of SU(u) Lie algebra basis is classified by total spin. In particular, for 
the S£/(3) problems, the primary with the identity matrix (A=0) is spin-singlet, three pri- 
maries with the antisymmetric matrices are spin-triplet and the remainder with symmetric 
ones are spin-quintuplet. The universal coefficients characterizing the S?7(3)-symmetric and 



asymmetric model are shown in Table III 



Let us now consider the third case g 2 < 0, which corresponds to 6 < 6 C = 7r/4 in the 
S = 1 model. The theory is asymptotically free, then we expect the mass generation which 
can be identified with the Haldane gap in the S = 1 case. One can estimate the mass gap 
by solving the renormalization group eq( |3.3|) . The conservation law eq (|3.4|) enables us to 
reduce the simultaneouse equation for the two unknown functions g\ and g% to that for the 
one unknown Y = —g\ — g% 



Ay / 

— = ±uyM^lY 2 VI + CY-^+ 2 )l\ (3.10) 

(JjI 

where the sign of the right hand side is identical to that of X = g\ — g^- Let us set the 
initial condition of the running coupling constants near the transition point 

<7i(0) = -oi </ 2 (0) ~a 2 {6-e c ), 

where a% and a 2 are positive constants. This condition sets the integral constant as C ~ 
a 3 (9—6 c ) in eq (|3^) with a positive constant 03. The renormarization group equation eq( p.!0|) 
is immediately integrated under this condition 
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where a = This gives us the order of the scale m 1 which makes the running coupling 
constant diverge g 2 {hi'm~ 1 ) = oo This scale m is the energy gap 

m = exp (-A\C\- a ) ~ exp (-c(6 c - 6)-°), (3.12) 

where 

a = 2 r°° d y 

V\JV 2 - 1 h y 2^Jl _ y-l/a' 

and c = a 3 A is positive. Therefore we conclude that the phase transition is infinite order. 
This result agrees with the recent numerical studies of S — 1 model by Fath and Solyom 
fl8fl . To see this, one should check their obtained energy gap directly rather than one 
parameter beta function estimated from it, since we have two parameter beta function fl3.3| ). 
Their numerical data of the energy gap fit the function eq( [3.12j ) with the universal constant 
a = 0.8 ± 0.2. This is consistent with our result a = = 0.6 at v = 3. 

IV. DISCUSSION AND OPEN PROBLEMS 

We have investigated the isotropic spin-1 model to clarify the phase diagram around the 
Uimin-Lai-Sutherland (ULS) point. The low energy theory of the ULS model is described 
by a strong coupling abelian gauge theory which can be regarded as the critical level-one 
SU(3) WZW model. 

We have shown a mechanism of the dynamical mass generation in the 5 = 1 Haldane 
phase in the presence of the SZ7(3) breaking interaction with dimension 2. We have shown 
that the dimension 2 operator makes the massless phase 9 < ir/A and the massive phase 
9 < 7r/4 around the ULS point 9 = tt/4 in the model eq ( PL . 1|) . This nature can be under- 
stood by the level k — 1 WZW theory, which has neither relevant operator with the chiral 
Z 3 invariance nor tensored operator of the WZ matrices but merely marginal operators. 
Therefore, the Haldane phase has the exponential mass gap as a result of the BKT tran- 
sition. The region 7r/2 < 9 < 7r/4 is concluded to be massless from this analysis and the 
numerical study |TB[ . Here, we indicate the difference of the phase transitions at the ULS 
point and at another integrable point 9 = — 7r/4 of the Takhatajan-Babujian (TB) model. 
In alternative field theoretical approach for understanding the Haldane massive phase, Af- 
fleck and Haldane investigated the relevant deformation of the S = 1 T-B model ||35|| . The 
universality class of this T-B model is the level-two SU(2) WZW model, where the one-site 
translation corresponds to the chiral Z 2 transformation. In the level- k theory with k > 1, 
one can make the chiral Z 2 invariant relevant operator in terms of tensoring of the SU(2) WZ 
matrices G(z,z), for example (Tr[G]) 2 . Therefore the transition from that massless point to 
the Haldane phase becomes second order, and the mass gap opens obeying the power law. 
In this case, the T-B point 9 = —tt/4 is isolated as a massless point in the massive region, 
namely the Haldane phase 9 > —tt/4 and the dimer phase 9 < —it /A. 

The renormalization group flow given by eq( ^.3| ) has a unique fixed point in v > 2 case, 
while that in v = 2 case has a fixed line. Contrary to the v = 2 case, the logarithmic correc- 
tions appears in the massless phase for v > 2 even if there is the S£7(V) symmetry breaking 
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interaction. We have calculated coefficients of logarithmic corrections to the energies of the 
ground state and some excited states both in SU(u) symmetric and asymmetric models. We 
find the different coefficients in these two cases. from their numerical data of the energy gap 
as in the form eq (|3.12j ) This nature of the model with v > 2 suggests Cardy's argument that 



a natural irreducible CFT with one parameter should have the central charge c vir = 1 []36 | . 

Nonetheless, no one has ever succeeded in classifying CFT with c vir > 1, and therefore to 
search CFT with fixed line (or surface) might be worth attempting. Since we need to spread 
the coupling constant space at least, the simplest candidate is a model with anisotropic 
parameters or q-deformation of the Lie algebra S[/(V). This program is now in progress. 

Here we present some conjectures deduced from the CFT kinematics. We note that J L (z) 
and J R (z) which are in a subalgebra of SZ7(3)i Kac-Moody algebra except the normalization 
satisfy the level-four SZ7 (2) Kac-Moody algebra. The representation of SZ7" (3) 1 is involved 
in that of S£/(2) 4 . The central charge of both theories are c vir = 2 and conformal 

weight of the primary field with spin-j is = j(j + l)/6 with < j < 2 [31]. If we 



neglect primaries with half-odd-integer spin in the S?7(2)4 WZW model, we obtain those in 
the S?7(3)i WZW model. The SZ7 (2) 4 WZW model can be regarded as a critical theory of 
the spin-2 TB model, and therefore we can expect the following prediction: 

Conjecture 1 There is a cross-over flow from the spin-2 Takhtajan-Babujian model to the 
spin-1 Uimin-Lai-Sutherland model. 

As recognized in the studies of the S£/(2) spin chains, coefficients 

dj in the logarithmic correction to the excited states with total spin-j satisfy the following 
sum rule |3^| : 3d t + ld s = 0, where d s ,^ is the universal coefficient for the singlet (triplet) 

excitation(s) and the prefactor is the dimension of the spin representation. We have seen 
that such a similar rule exists in the spin-1 models discussed 



above, as well. That is 5d + 3d t + d s = f3"3] . Therefore, we are led to the following 
conjecture: 

Conjecture 2 There exists a sum rule among the coefficients {dj} of the leading logarithmic 
correction term in the excited energy with total spin-j ; i.e. 

2S 

]T(2j + i)^ = 0. 
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APPENDIX A: 



The fundamental representation of the SU(u) Lie algebra [T , T ] = f C T is sum- 
marized as follows. The Sf/(z/) exchange operator is decomposed in terms of SU(u) basis 
as 

1 V2 - X 

V = -I x / + 2 Y T A x T A . 
These basis are normalized as Tr [T A T B ] = ^5 AB 



or 



Ea=i T a T a 



l)/2v. The structure constant f ABC has the quadratic Casimir 



of the adjoint representation: J2a~b=i f ABC Iabd = ~ v ^° d- Another expression of the 
exchange operator is available when the spin chains are studied. On a space C 2S+1 x C 2S+1 , 
it is given by 



V 



as 



25 



E(-i) 3 ^. 



3=0 



where V^' is the projector onto a space of spin-i conforming to an identity 1x1 = + 
• • ■ + T >{ ' 2S \ The projector on a spin-j space is represented using the spin operators with 
the magnitude S as follows: 



V 



2S 

n 

fc=0 
&3) 



X - x k 



X i 



X = s a x S a , 



a=l 



where x k = [k(k + 1) — 2S(S + l)]/2. The expressions of the exchange operator in terms of 
the spin operator are shown in Table |. 

In particular, the representation of SU (3) is realized by Gell-Mann matrices 




where T A = X A /2. Here A j4=2 5 7 are antisymmetric matrices and the remaiders of them are 
symmetric. 

The primary states {\O a )} can be classified by total spin-j. The total spin operator is 
given by 

^tot — dx S(x) = J L + J R0 , 
Jo 



where SZ7 (2) charge operators are J L0 



~J L {z) and J R:0 



■ dz 
2yri 



J R {z). The magnitude of 



total spin of the primary states takes values 0, 1 or 2 from a synthesis of two fermions with 
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spin I. Acting S tot to the primary fields, we obtain the OPE 
S tot O A =°(z,z) = 0, 

S tot O A ^(z,z) = AO A (z,z) + 2T A M*(z)^(z)e i ^ ) + h.c). 

Here we have used the properties of the SU (3) basis. Using the symmetric and asymmetric 
properties of the Gell-Mann matrices, we obtain 

(s tot y\o A )=j(j + i)\o A ), 

where j = 0,1,2. The primary with identity matrix (A = 0) is the singlet state (j = 0). 
Three antisymmetric ones (A=2,5,7) in the Gell-Mann matrices give spin-triplet states (j = 
1). The remainders which are symmetric matrices, become quintuplet states (j = 2). 
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TABLES 



TABLE I. Expressions of the exchange operator for S < 2 



s 


V 




1 

2 

1 

3 
2 

2 


2X + i 
X 2 + X - 1 

2 v3 i 11 v2 
9 A + T8 A ~ 
1 Y 4 _i_ 1 ^3 
36 + 6 A ~ 


9v 67 
8^ 32 
7 Y 2 5 Y 1 
T2 A ~ 2 A ~ 1 



TABLE II. Finite-size corrections for the spin-1/2 Heisenberg chain 



c vir 


x t x s d GS 


dt 




d s 




SU(2) 1 WZW 1 


1/2 1/2 3/8 


-1/4 




3/4 


1/2 


BA 1 


1/2 1/2 0.3433 


-1/4 




3/4 




TABLE III. 


Finite-size corrections for the 


spin-1 chains 








Wir 


x q x t x s d GS 




d t 


d s 


°? 


SC/(3)i WZW (7 = 1) 2 


2/3 2/3 2/3 4/9 


-1/9 


-1/9 


8/9 


2/9 


SC/(3)i WZW (7/I) 2 


2/3 2/3 2/3 6 


-1 


1 


2 


2 


BA 2 


2/3 2/3 2/3 
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FIGURES 



FIG. 1. The renormalization group trajectory for the v > 2 model, where v = 2S + 1. The 
coupling g 2 is defined in the vicinity of the ULS model. The BKT line g 2 = 0, g 1 < corresponds 
to the pure ULS model. 
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